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Abstract 

We study the complexification of the one-dimensional Newtonian 
particle in a monomial potential. We discuss two classes of motions on 
the associated Riemann surface: the rectilinear and the cyclic motions, 
corresponding to two different classes of real and autonomous Newto- 
nian dynamics in the plane. The rectilinear motion has been studied 
in a number of papers, while the cyclic motion is much less under- 
stood. For small data, the cyclic time trajectories lead to isochronous 
dynamics. For bigger data the situation is quite complicated; com- 
puter experiments show that, for sufficiently small degree of the mono- 
mial, the motion is generically periodic with integer period, which de- 
pends in a quite sensitive way on the initial data. If the degree of the 
monomial is sufficiently high, computer experiments show essentially 
chaotic behaviour. We suggest a possible theoretical explanation of 
these different behaviours. We also introduce a one-parameter family 
of 2-dimensional mappings, describing the motion of the center of the 
circle, as a convenient representation of the cyclic dynamics; we call 
such mapping the center map. Computer experiments for the center 
map show a typical multi-fractal behaviour with periodicity islands. 
Therefore the above complexification procedure generates dynamics 
amenable to analytic treatment and possessing a high degree of com- 
plexity. 
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1 From Newtonian systems to the motion on 
Riemann surfaces 

In this paper we investigate the following two classes of autonomous Newto- 
nian dynamical systems in the plane 



x = F^ m) (x), m E Z, ( | , 

x = F^^x, x), m G Z, m/-l, 



where x = (x,y) E M 2 , and / = d 2 f/dt 2 . The forces Fj ™ are defined by 
Fi m) (x) = - (Re(a 2 (x + iy) m ), Im(a 2 (x + iy) m ))) , 

F< m) (x, x) = —x A h m — 87r 2 ^±^x - (Re(x + iy) m , Im{x + iy) m ) , (2) 

where the constant vectors h m , m G Z are orthogonal to the (x, y) plane, 
with ||hj| = 2vr^±f, and a G C is an arbitrary parameter. 

The classes (J^i) and (JTJ)) originate from the same complex ODE: 

= -C mGZ, m^-l, TEC, (e((t)6C, (3) 

solved by the complex quadrature 

C 

C( J 0) V^¥)' (4) 
^=Hf(°)) 2 + ^(C(0)r, n = m + l. 

The first class ((Ui) is directly obtained from (jHJ) via the following change of 
independent variables 

z(t)=x{t)+iy(t) = ((r), , , 

r(t) = at, a EC { ' 

corresponding to a rectilinear motion on the hyperelliptic curve of the quadra- 
ture dU); the second class (JTJd) is instead connected to (jSJ) via the following 
change of dependent and independent variables (Calogero's transformation) 

DP: 

z (t) = x(t) + iy(t) = exp (i±g) C (r) , 

/,S __ Cxp(t2 7Tf)-l l D J 

v ' 2ni 



Equation (J^b) implies that, as the "physical time" t travels onward along 
the real i-axis, the motion of the complex time r is cyclic; the auxiliary 
transformation (jHK) is what makes the Newtonian class (03 ) autonomous. 
Through the change of variables 

w(0 = c iC(t), £ = c 2 t + £ , 

ciC(o) (7) 



dw 







the quadrature (J3J) takes its adimensional form. We have 2 cases, correspond- 
ing to n positive and n negative respectively: 



u, (8) 



dw / w\ n \ 
oj = —== , n > 0: oj = \ — rn dw, n < 0. (9) 

vi -w n v wH - 1 

Therefore the Newtonian systems (JTJi) and ([TJd) correspond, via (J3j), (0), 
(0), respectively to the rectilinear motion in the £-plane: 

£ = £ + «c 2 t (10) 
and to the cyclic motion in the £-plane: 

£ = £ c + #exp(i(27rt + ^)), 

e c = e -^ ? ^ = arg(-zc 2 ), 1 j 

where £ c and R are the center and the radius of the corresponding circle. 
We remark that the Cauchy data of the Newtonian dynamics (0d) determine 
completely this circle and the initial point on it. 

To express w as a function of t we have to invert the hyperelliptic integrals 
©, © defined on the surfaces 

r : fi 2 = 1 - w H , n > or n = 2k < 0, (12) 
r x : /i 2 = w - w H+1 , n = 2k + l<0; (13) 
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therefore the trivial dynamics (|1U|) and in the £-plane result in highly 
nontrivial motions on T, Fx. 

Let us point out that, for n > 0, we have the problem of inverting a 
hyperelliptic integral defined by a holomorphic 1-form, and for n < we 
have to invert an integral defined by a meromorphic 1-form. As we shall 
see below, if the motion is rectilinear, the case n < is much simpler from 
the dynamical point of view. 

The rectilinear motion on Riemann surfaces can be naturally interpreted 
as Hamiltonian system with a multivalued Hamiltonian ( or, equivalently, 
as a Hamiltonian system defined by a Hamiltonian 1-form with non-zero 
periods, see review (THj). This motion was intensively studied in the literature 
(references can be found in the reviews |24[ II lj). A new description of the 
topology of these flows was suggested recently in [T5] . 

The interest to this problem was partially motivated by its connections 
to stochastic dynamics. Consider a cycle on the Riemann surface transversal 
to the flow. The Poincare map on this cycle preserves the measure on the 
cycle, therefore it can be interpreted as the interval exchange map, which 
is one of the basic stochastic models. Moreover, any interval exchange map 
can be obtained from this construction. It turns out that the properties 
of this map are deeply connected with the geometry of the moduli spaces 
O H3 122 I2H EH]- In particular, the Lyapunov exponents for this map 
at generic points can be calculated in terms of the Lyapunov exponents for 
the geodesic flow on the Teichmuller space. 

Another motivation for the study of rectilinear motions comes from solid 
state physics: the conductivity tensor for metals in a "strong enough" mag- 
netic fields is deeply connected to the quasiclassicl electron motion on the 
Fermi-surface, which is a Hamiltonian flow on the Fermi-surface, and the 
Hamiltonian 1-form is determined by the direction of the magnetic field > 
jllj . In contrast to the case of the interval exchange map, the holomorphic 
form is now non-generic, and its real part has exactly 3 non-zero periods; 
therefore this situation needs a separate study. 

The study of cyclic motions on Riemann surfaces was initiated in pQ, 
where a general procedure, involving the change of variable (jE|L was intro- 
duced to deform large classes of ODEs so that they feature a lot of periodic 
solutions. The deformed dynamical systems are isochronous for initial data 
lying in a certain domain of full dimensionality (corresponding to the situ- 
ation in which the circle is contained inside the analyticity disc of the un- 
modified ODE). Outside this domain, the circle may include branch points; 
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then bifurcations occur, and either the motion becomes periodic with mul- 
tiple period, or the periodicity is lost and the system may exhibit sensitive 
dependence on the initial data (see 0121013111131110 and references 
therein quoted). This mathematical mechanism was first described in j3Cj- 
As it was shown in [S] for an aristotelian 3-body problem which allows for 
a quite explicit analytical treatment, the transition from regular (periodic) 
to irregular (with sensitive dependence on the initial data) dynamics takes 
place when an infinity of branch points of the solution are arbitrarily close 
to the real t-axis. 

A very general class of autonomous dynamical systems obtained through 
Calogero's transformation has been constructed in |S] (see also |15j). This 
class contains, in particular, systems of ODEs solved by quadratures asso- 
ciated with elliptic and hyperelliptic curves; the first investigation of some 
examples of this type was presented in [T3j . Also the Newtonian systems jib), 
corresponding to cyclic motions on the Riemann surface of the quadratures 
flU), are examples belonging to this general class. A detailed study, analytical 
and computational, of the systems (JUd) in the particular cases n = 4, 5, 6 has 
been recently done in ^2], with the following results. In the elliptic case 
n = 4 the motions are isochronous with period 1; if n = 6 the hyperelliptic 
curve doubly covers an elliptic curve, and all motions are periodic with peri- 
ods 1 or 2. In the more interesting case n = 5, numerical experiments seem 
to indicate that all solutions are periodic, with arbitrarily large period. The 
results contained in (3 El H2] have been the main motivation of our study. 

Remark 1. The cyclic motion can be described as the motion with con- 
stant geodesic curvature on the Riemann surface with almost everywhere 
flat singular metric (or, equivalently, as the motion with a constant geodesic 
curvature on a polyhedron). As it was pointed out by S.P.Novikov ^H], i n 
the 2-dimensional case the motion with a prescribed geodesic curvature can 
be naturally interpreted as the Hamiltonian motion in magnetic field. 

Remark 2. Under the change of variables (J3), the complex Newtonian equa- 
tion © becomes z = —a 2 z m , and can be cast into the complex Hamiltonian 



Then it follows that also the corresponding dynamical systems (JUi) can be 



form: 




C 



(14) 
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written in Hamiltonian form: 



i — — 

dpx ' 




(15) 




dH 

dy ' 



where the canonically conjugated real variables x = (x, y) and p = (p x , p y ) are 
defined by z = x + iy, it = p x — ip y , and the Hamiltonian H(x,p) is the real 
part of 7i (the general procedure to go from complex to real Hamiltonian 
formulations can be found in [21] )■ In addition, the real and imaginary 
parts of TC are constants of motion in involution for the Hamiltonian system 
()15|). Therefore the Newtonian equations (^i) are Liouville integrable, with 
polynomial (if n > 0) or rational (if n < 0) constants of motion in involution. 
But this Liouville integrability is not global. Indeed, a simple local analysis 
shows that the singularities of x(t) and p(t) in the complex t - plane are, for 
n 7^ 3,4, branchpoints (see (J2U|0 . In addition, for n^6, such branch points 
are everywhere dense in the complex t - plane, x, p are unbounded and 
the 2-dimensional real variety defined by the two constants of motion is not 
a torus (in particular, global action - angle variables cannot be defined). 
The existence of everywhere dense branch point singularities, which causes 
the lack of global integrability in the Liouville sense, turns out to be the 
reason for the sensitive dependence on the initial conditions exhibited by the 
corresponding dynamics, discussed in the following sections. 

For the Newtonian class (0d), which clearly inherits from Q and © the 
two t-dependent first integrals H, K: 



no Hamiltonian structure is instead known. Only after fixing the energy 
level, one obtains the Hamiltonian interpretation contained in the Remark 
1. 

Remark 3. The idea that the integrability or the non- integrability of a dy- 
namical system be intimately related to the analyticity properties of its so- 
lutions in the complex time plane goes back to Jacobi, Kowalewski, Painleve 
and other eminent mathematicians. For Kowalewski, Painleve and his school, 
the request that the only movable singularities of ODEs be poles has been 
a very successfull tool to isolate and classify important integrable cases ^3] . 
The idea that integrability is compatible even with the presence of movable 



H = Re H, K = Im H, 



(16) 



12^ n— 2 / n J 



- 47rn 
n-2 
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branch point singularities, provided that they are not dense in the complex 
time plane, is more recent, and one of us (PMS) heard several illuminating 
lectures on this topic by Kruskal, more than 20 years ago. The works 120] 
are some of the contributions connected to this fascinating idea. 

As we shall see in the following sections, the Newtonian systems (JQ are 
just distinguished examples of ODEs whose solutions exhibit branch point 
singularities which are everywhere dense in the complex time plane. There- 
fore the results presented in this paper can also be viewed as a contribution 
along this line of thinking. 

The paper is organized as follows. In §2 we discuss the topological prop- 
erties of the Riemann surfaces associated with the quadratures for a 
generic n G Z. In §3 we concentrate on the physically relevant subcase in 
which the degree n of the monomial is even, restricting further our analysis to 
a suitably factorized Riemann surface. In §4 we discuss the main properties 
of the rectilinear and cyclic motions on such factorized surfaces, as the degree 
n of the monomial varies. In §5 we associate with the cyclic motion the 1- 
parameter family of 2-dimensional mappings describing the dynamics of the 
center of the circle; such mappings show typical multi-fractal behavior with 
periodicity islands. In §6 we briefly comment on how the complexification 
procedures discussed in this paper can be a distinguished way of generat- 
ing dynamics exhibiting a high degree of complexity and, at the same time, 
amenable in principle to analytical treatments. 

2 Riemann surfaces associated with the New- 
tonian dynamics: generic case 

In this section we introduce a description of the Riemann surfaces associated 
with the quadratures ©-©• This description is different from that used in 

m 

n > 0. In the case n > we use the following representation for the surface 
T. Due to the rotation symmetry of T, it is convenient to draw the cuts in 
the w-plane as rays connecting the branch points (the n roots of unity) to 
infinity. Then we have 2 sheets and the map (jHJ) maps each of them to a 
regular n-gone (see Fig 1): 
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w-plane 



\ - plane 




Figl. 

Therefore the whole Riemannn surface is represented by the union of 2 
regular n-gones, and each side of the first n-gone is glued to the opposite 
side of the second one. Equivalently, we can attach one side of the first n- 
gone to the corresponding side of the second one. As the result of the above 
operation, we obtain a {2n — 2)-gone, and Y is obtained by gluing the opposite 
sides (see Fig 2). 

(2n-2)-gone, n=8 




Fig 2. 

Let g(T) be the genus of T; then, for both n 
g(V) — k — 1. The length of each side is equal to 



2k and n — 2k — 1, 



I 



dw 



\/w n - 1 



(17) 



8 



and the distance h from the center to the sides satisfies the relation I = 
2/itan(7r/n). 

If n is even, it is also convenient to define the factorized Riemann surface 
T, obtained by the factorization 

w — > — w, /i — > — fx. (18) 

T can be naturally treated as the result of gluing the opposite sides of a single 
regular n-gone ^H]- H n = 4k or n = 4k + 2, then g(T) = k. 

n < 0. In the case n < it is convenient to draw the cuts as rays connecting 
the roots of unity to 0. If n is even: n = 2k, we have 2 sheets and the 
map Q maps each of them to the exterior (complement) of the regular \n\- 
gone. Therefore the whole Riemannn surface T is represented by the union 
of the complements to 2 regular |n|-gones, and each side of the first |n|-gone 
is glued to the opposite side of the second one. Topologically, the exterior 
of an n-gone with the infinite point added is an n-gone; therefore T can be 
represented again as the result of gluing the opposite sides of a (2n — 2)-gone. 

If n is odd: n — 2k + 1, we have one sheet and the map © maps it to 
the exterior of a regular (2|n|)-gone, which is drawn on the Riemann surface 
of the function ^Jw. The points of a side in one sheet are glued to the 
corresponding points of the side in the second sheet. 

Cycles. In our situation it is more convenient to use, instead of the canonical 
basis of cycles, a more symmetric one. Since V is obtained by gluing the 
opposite sides of a (2n — 2)-gone, let us choose the center of this {2n — 2)- 
gone as the starting point and denote by Cj the cycles connecting the sides 
Sj and s J+n _i, j = 1, . . . , n — 1 (see Fig 3). 
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(2n-2)-gone, n=8 




If (n — 1) is even, these cycles are homologically independent, Cj • Cj — 1 
for i < j. The fundamental group of T is generated by the Cj with the 
relation dC^C^ 1 . . . C^ 1 C^ 1 C 2 C^ 1 C 4 . . . C n _ x = 1. 

If (n — 1) is odd, these cycles are dependent. It is convenient to write 
the fundamental group of T as the group generated by Cj, j — 1, . . . , k with 
two relations: dC^C^ 1 . . . C n _i = 1, C^dC^Ci . . . C"^ = 1. In the 
homological group these two relations are equivalent and they imply that 

C 1 -C 2 + C 3 -C 4 ... + C n ^ = 0. (19) 

As usual, the multivaluedness of £(w), associated with the above cy- 
cles, implies that the inverse function w(£) be mult iper iodic. Moreover, the 
vertices of the n-gone are the only singularities of w(C,). If n = 3,4, this 
singularities are poles and is single- valued. Otherwise, they are branch 
points £ 6 , and w(£) is unbounded near them. For n > 

w~c(Z-Q-^, (20) 

Due to the above multiperiodicity, these singularities are everywhere 
dense in the complex plane for n ^ 3,4,6. These exceptional cases corre- 
spond to the equilateral triangle, the square and the regular hexagon, which 
are well-known to cover exactly the plane through the Schwarz reflection 
principle. 

Coverings. The inversion of the hyperelliptic integrals leads us to the study 
of the motion not only on the surface T, but also on a proper covering space 
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f . It is well-known that there is a one-to-one correspondence between the 
subgoups 7f of the fundamental group n(T, x ) and the non-ramified coverings 
of T. Namely, the paths 7 G tx are exactly the loops which remain closed 
after being lifted to covering space. If % is the unit subgroup, we obtain the 
universal covering, which for g > 2 is isometric to the Lobachevsky plane. 
But, for the problem of inverting the hyperelliptic integrals, it is much more 
natural to use "smaller" coverings: 

1. The maximal Abelian covering (see JH|): 7r is the commutant of tt, 
i.e. the subgroup generated by all elements of the type g\g29\ X g^ 1 > 
9ii92 £ 7r(r,:co). A closed path in Y remains closed in the covering 
space iff the integrals of all holomorphic forms (as well as the integrals 
of all meromorphic forms with only zero residues) along this path are 
equal to 0. 

2. The covering T associated with a given holomorphic form u: the sub- 
group 7r is formed by all closed paths 7 such that 

j u = 0. 

7 

If all periods of u are linearly independent over the rational numbers, the 
covering defined by u coincides with the maximal Abelian one. But, if there 
is some rational dependence between the periods, it is smaller. 

3 Riemann surfaces associated with the New- 
tonian dynamics: n even 

For simplicity, in the rest of the paper we concentrate only on the case of 
even n: n = 2k. We shall also impose the restriction n > 0. The reason is 
that for n < the rectilinear motion is trivial (see below), while the cyclic 
motion can be obtained from the case n > by time- reversal. 

It is convenient to focus our attention on the motion on the factorized 
Riemann surface T. If the motion on T is periodic, then the motion on F 
is also periodic, and the two periods either coincide or the period on T is 
twice bigger, depending on whether the number of intersections with the 
sides is even or odd. Therefore the motion on T contains all the essential 
informations. The corresponding cycles are drawn in Fig 4. 
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Fig 4. 



Since our experiments cover the polygons with even n, from n = 8 to 
n = 16, let us discuss the main properties of the corresponding Riemann 
surfaces F, F. 

If n = 4, T and F are tori: g(F) = g(F) = 1, and the universal covering 
coincides with the maximal Abelian covering and is exactly the complex 
plane. The basic periods are 2h and 2ih, and they are independent. 

If n — 6, f is again a torus, g(F) = 1, g(F) = 2, and T is a ramified 
covering of T with 2 branch points. The periods along the cycles C±, C 2 , 
C3 are 2h, 2e tn ^ 3 h, 2e 2t7T ^ 3 h respectively, and their alternated sum is equal 
to 0. Let us recall that, homologically, C 1 — C 2 + C 3 = 0. The covering T 
corresponding to our integral coincides with the maximal Abelian covering. 

If n — 8, g(T) = 2, g(T) = 3 and T is a non-ramified covering of T. The 
basic periods of T are 2h, 2e l7T ^ 4 h, 2e tn ^ 2 h, 2e 3ln ^ 4 h, and they are independent 
over the rational numbers. Therefore the covering T associated with our 
integral again coincides with the maximal Abelian covering. 

If n = 10, g(T) = 2, g(T) = 4 and T is a ramified covering of T with 
2 branch points. The basic periods of f are 2h, 2e l7r/5 h, 2e 2i7r/5 h, 2e 3n/5 h, 
2e 4m,/5 /i, and the basic relation is C\ — C 2 + C 3 — C 5 + C 5 = 0. Again the 
covering F associated with our integral coincides with the maximal Abelian 
covering. 

If n = 12, g(F) = 3, g(F) = 5 and T is a non-ramified covering of f\ 
The basic periods of f are 2h, 2e iw/6 h, 2e iw/3 h, 2e iw/2 h, 2e 2iw/3 h, 2e 5iw/6 h; 
therefore 

j cu = j cu = 0. (21) 

Ci — C3+C5 C2 — C4+CQ 
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It is the first case in which we have nontrivial integer relations between the 
periods: we have 6 independent directions in the universal covering, but 
only 4 of them remain independent on the covering T associated with our 
hyperelliptic integral. 

If n = 14, g(T) = 3, g(T) = 6 and T is a ramified covering of T with 2 
branch points. The only relation between the periods is the one prescribed 
by the homological group: C\ — C 2 + C 3 — C 4 + C 5 = 0. The covering T 
associated with our integral coincides with the maximal Abelian covering. 

The last case studied in our numerical experiments is n = 16. In this case 
g(T) = 4, g(T) = 7 and T is a non-ramified covering of T. The basic periods 
of f are independent over the rational numbers, and we have the full Abel 
covering. 

The points of the maximal Abelian covering are represented by the points 
of the basic polygon plus the total shift 

k 

T = J2m J C J , (22) 

3=1 

where rrij is the number of crossings of the trajectory with the side Sj minus 
the number of crossings with the opposite side Sj + k- For odd k, the relation 
(fTU|) should be taken into account. The simplest way is to substitute Ck = 
-C\ + C 2 + ... + C k -i into (J22), obtaining 

fc-i 

T = EK' + ( 23 ) 

3=1 

Due to the rotation symmetry of T, we can have extra integer relations 
among the periods of the form. For instance, if n = 41 k' , where I is odd, 
the cycles C m , C m+ 2k', Cm+ik', ■ ■ ■Cm+2k'(i-X) are independent, but the cor- 
responding periods satisfy the relations: 

cj = 0, m= 1,2..., 2k'; (24) 

see, for instance, (j2"Tj) for n = 12. In these cases, instead of studying the 
maximal Abelian covering, we study the smaller covering T associated with 
our hyperelliptic integral, and we have extra relations. Therefore 

h 

T = Y,rhfi 3 (25) 
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for suitable rhj, where n is the number of independent periods. 

4 Rectilinear and cyclic motions 

Rectilinear motion. If n < the qualitative picture is rather simple. The 
integrand is a meromorphic differential, and each trajectory goes to oo for 
t — > ±00. In terms of the vector field, it means that we have either 2 (for 
even n) or 1 (for odd n) fixed points, all the trajectories reach them in an 
infinite time, and generic initial conditions result in regular motions without 
sensitive dependence on the initial condition. 

If n > we are in the situation in which the Hamiltonian flow is defined 
by a holomorphic form. The basic features of such motion for g > 1 were 
intensively studied by many authors (see the literature at the end of §1). It 
is convenient to choose a Poincare cycle transversal to the flow, and study 
the corresponding Poincare map. For a given trajectory, the number N of 
intersections with the Poincare cycle is used as the discrete time. 

For generic data (a generic algebraic Riemann surface and a generic holo- 
morphic form) the motion is aperiodic, the image on the maximal Abelian 
covering has a distinguished direction v (which is Poincare dual to the co- 
homological class of the form). In our case 

k 

vo = Y^vdjCj (26) 

where 

(vo)j = (a,fij), a = (Rea,lma),hj = (cos(27r(j — l)/n, sin(27r(j — l)/n)), 

(27) 

■hj denotes the normal vector to the side Sj. Let us recall that, if n — 2k 
and k is odd, the vectors Cj are dependent, and for the maximal Abelian 
covering we can use the basis Ci, . . . , Ck-i- Then 

fc-i 

^o = $>o),Q, («b)i = Wi + (28) 
3=1 

If the number iV of iterations is sufficiently large, the total shift from the 
starting point is described by the formula 

T = Ncv + o(N), (29) 
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where c is an order 1 constant. The Poincare map is ergodic, and this implies 
that, for any pair of arbitrary close initial points, the stripe between their 
trajectories sooner or later meets a branch point, located at the vertex of the 
polygon. After this event the trajectories separate to a distance of order 1. 
The difference between these trajectories is encoded in the o(N) term of the 
asymptotics. The corrections to the leading asymptotic term were studied 
in j2H|, p6j (see also the review [21]), where the following statements were 
proved: 

1. Let us consider the projection P^ along the vector v to the comple- 
mentary space. Then there exists a vector v\ in this complementary 
space such that the projection has the form: 

p^T = N Xl Cl (N)vi + o(N Xl ), < Ai < 1. (30) 

where C\ (N) is a random function of N of order 1 

2. For any given starting point, it is possible to define a special sub- 
sequence {N t } such that, going from Ni to N L+1 , the trajectory re- 
turns much closer to the starting point (for detailes see fH L ). Using 
this subsequence, it is possible to continue the above procedure, defin- 
ing the vectors v^, ■ ■ ■ , v g -\ and the sequence A2, . . . , A 9 _i such that 
1 > Ai > A 2 ...A fl _i > 0. 

3. If o denotes the intersection form on cycles, it generates a symplectic 
scalar product on vectors. In this symplectic product Vi o Vj — for all 

If we consider a pair of close trajectories, after a sufficiently large number 
of iterations they deviate on the maximal Abelian covering, and, according to 
the formula (J3U|) . the distance on the maximal Abelian covering increases as 
iV Al . Therefore Ai can be interpreted as the main Lyapunov exponent, and 
the remaining A/s can be interpreted as the other Lyapunov exponents (see 
24J). Additional information about the numerical values of these exponents 
for small g as well as some analytic results can be found in the review |24j . 

Formally, the results of the papers [2I]-i2II| are valid for generic Riemann 
surfaces, while our polygon is highly symmetric. But numerical experiments 
show good agreement with this theory. 
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Cyclic motion. As far as we know, the study of cyclic motions on Riemann 
surfaces began only recently (see the literature at the end of §1) and the 
picture is much less clear. Here we present the results of numerical experi- 
ments proceeded by the authors, and we give some preliminary qualitative 
explanation of the observed phenomena, based on the topological structure 
of the surface. 

We always assume that our trajectories are generic in the following sense: 
they never hit the vertices of the polygon, otherwise the motion becomes 
undefined. 

For the numerical simulations we used the following technique. As it 
was mentioned above, we make our calculations on T, which is represented 
by gluing the sides of the regular n-gone, n = 2k. In fact, to recover the 
dynamics on T for even n, it is sufficient to study the dynamics on T and 
check at each step if the number of intersections of our trajectory with the 
sides of the polygon is even or odd. If it is even, the trajectory is on the 
same sheet, and if it is odd, the trajectory is on the opposite sheet. 

Any time our trajectory crosses the side of the polygon, it is shifted of 
the corresponding period. 

In the cyclic case the trajectory is periodic iff, at some iteration, the total 
shift on the covering associated with the quadrature becomes zero. 

We first discuss the exceptional cases. If n = 4, T and T are tori, w(£) 
is single-valued and the motion is always periodic with period 1; i.e. it is 
isochronous. If n = 6, T is a two-sheeted ramified covering of the torus T with 
two branch points. Therefore the motion is always periodic, with period 1 or 
2, depending on whether the number of square root type branch points inside 
the trajectory is even or odd. These two cases have been already considered 
in [ig. 

We now discuss the non-exceptional cases. 

For given n = 2k, if the diameter D of the circle is smaller than /, where 
I is the side length, the motion is periodic, and there is at most 1 branch 
point inside the circle. If there is no branch point inside, the period is 1 and 
the motion is isochronous; otherwise it is equal to {n — 2)/2 for even k and 
to (n - 2)/4 for odd k. 

For bigger cycles the situation becomes much more complicated, and we 
studied it using a series of computer experiments, which have not exhausted 
the whole phase space. These experiments show the following. 

1. The period is very sensitive to the radius R and to the position £ c of 
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the center of the circle. 

2. Our observation shows an essential difference between the cases n = 
8, 10, 12 and n = 14, 16. 

Let us discuss the first group n = 8, 10, 12. From the numerical exper- 
iments it seems that, for data of full measure, the motion is periodic, 
and the typical periods grow exponentially with the ratio D/l. We 
have observed periods of order 10 8 for D/l ~ 20. For bigger D, we 
have observed starting points for which the period was not found after 
2 • 10 9 iterations. 

For n — 14 and n = 16 we have found that, for many initial configu- 
rations starting from D/l ~ 2, the period was not found after 2 ■ 10 9 
iterations. Our conjecture is that, in this cases, the aperiodic trajecto- 
ries are non-exceptional. 

3. If we consider the trajectories for which the period was either not found 
or was sufficiently large, the total shift on the covering space seems to 
exhibit the following behaviour: 

T ~ N x v(N), A ~ -, (31) 

where N is the large number of iterations and v(N) is a random vec- 
tor of order 1. Therefore, the distance on T of two close trajectories 
increases approximately as N x , and A ~ ~ can be interpeted as the 
principal Lyapunov exponent. 

Formula (|3T|) suggests that, for sufficiently large N, the shift T can be 
approximated by the random motion on a lattice. Now we show that, ac- 
cepting this assumption, one can explain the difference between the above 
two classes. 

Let us stress that, for cyclic motions, the total shift in the covering space 
can be arbitrarily large but, after being projected to the time-plane, the total 
shift is bounded by the condition: 

\J2 m j eU ~ 1)n/k \ < 4 + l/cos(7r/n). (32) 

This inequality imposes a severe constraint on the vector (mi, . . . , m&). 
Let us denote by "effective dimension" the maximal number of m/s which 
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can be chosen arbitrarily. If these components are fixed, there is only a finite 
number of options for the other components. 

Let us point out that the effective dimension is equal to the number of 
independent shifts in the covering F associated with the quadrature minus 2. 
Indeed, if we fix any pair of basic shifts, say e^' _1 ^ 7r /*, j = 1, 2, which form a 
basis in the complex plane over K, and we assign arbitrary integer values to 
all the other m/s (j ^ 1, 2), then the two coefficients m 1; m 2 are constrained 
by the inequality (|32|) to take only a finite number of values. 

From the considerations made in Setion 3, it follows that, for n = 8, 10, 12 
the effective dimension is 2, while for n — 14 it is equal to 4 and for n = 16 
it is equal to 6. 

On the basis of the numerical experiments we assume that, for sufficiently 
large N, our dynamics be approximated by a random motion in the "lattice 
of effective dimension". It is well-known that, for the 2-d symmetric random 
motion, the probability to return to the starting point is equal to 1; therefore 
we expect that, with probability 1, this motion be periodic. If the dimension 
is higher than 2, with non-zero probability we do not return to the starting 
point, and there are aperiodic trajectories. Therefore the random hypothesis 
implies that, for n = 8, 10, 12 (cases in which the effective dimension is 2) the 
typical motion is periodic, while for n = 14 and n = 16 (cases in which the 
effective dimension is greater than 2) there are aperiodic domains of positive 
measure, in agreement with our observations. 

5 The center map as a fractal 

In the rectilinear motion it is convenient to represent the dynamics on the 
surface in terms of the Poincare map on a transversal cycle. It is the well- 
known "interval exchange map". In the cyclic motion it is convenient to 
study the dynamics of the center £ c of the cycle, for any fixed radius R. 

The dynamical rule is the following: any time the circle crosses a side Sj 
of the n-gone, the center is shifted in the following way: 

£ c -> £ c - 2he 27T(j - 1)/n . (33) 

It is a 1-parameter family of deterministic 2-dimensional mappings (the ar- 
bitrary parameter being the radius of the circle). We call this mapping the 
"center map". 
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It is interesting to study the center map for aperiodic trajectories (of 
course, up to the computer experiment limits). Considering as basic illustra- 
tive example the polygon n = 16, h = 1 centered at the origin, we have made 
a series of experiments showing interesting generic multi-fractal behaviour 
with periodicity islands. 

In the first set of experiments we fixed R = 0.358 and we varied the initial 
position of the center. In Fig 5 we choose £ c = (1.1610, —0.4037), obtaining 
the following fractal: 




In this figure, together with the fractal, we see also the 16-gone and the 
initial circle, located in the south-east position (it is drawn by a very thin 
line). The initial center £ c is marked by the intersection of the inner lines. 
We remark that this fractal shows the rotation symmetry of the polygone. 

In Fig 6 we see the same fractal, but magnified 33 times. The initial 
center is marked again by the intersections of the inner lines. 




Fig 6. 



20 



This figure makes clear the fractal structure. It turns out that this fractal 
is the boundary of some islands inside which the center map is periodic 
with fixed periods. Typically the period increases as the size of the island 
decreases. 

Choosing the starting point £ c inside the fractal region, we obtain a fractal 
coinciding with this one up to the picture resolution. Choosing the starting 
point £ c = (1.1572, —0.4151) ouside the fractal, but on the boundary of the 
periodicity islands, we obtain a different fractal, living in the complementary 
region (see Fig 7). 




Fie 7. 
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In the second series of experiments we change the radius R. For a suffi- 
ciently small change of radius the fractal changes, but keeping all the prin- 
cipal features. Choosing, instead, a sufficiently different radius, we obtain 
completely different fractal pictures. For instance, for R = 0.6, we obtain 
the rich fractal shown in Fig 8 (compare it with the fractal in Fig 5). 



Polygons _ □ x 




Fig 8. 
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Fig 9 shows a detail of the fractal of Fig 8, magnified 10 times, containing 
a part of the initial circle. 




Fig 10 shows a detail of Fig 9, magnified 4 times more. These magnifica- 
tions make evident the extremely rich fractal nature of the center map. 




24 



6 Complexification and Complexity 



We end this work with a consideration. The results of this paper strongly 
suggest that proper complexifications of dependent and independent variables 
of a given ODE are very fruitful and rich operations, enabling one to generate 
complicated dynamics amenable to exact analytical treatments. Indeed they 
allow to go, for instance, from the elementary one dimensional Newtonian 
motion in a monomial potential, to the complicated dynamics (JU 3 ), whose 
high degree of complexity is well illustrated by the fractal behavior of the 
associated center map (see §5). We feel that dynamics generated through 
complexification procedures could serve as prototypical examples of a new 
paradigm of chaos, potentially interesting in applications. 
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